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Continuation of solutions of coupled dynamical systems^ 

Tianping Chen^, Wei WiO 
Abstract 



o 
o 

(N 

Recently, the synchronization of coupled dynamical systems has been widely stud- 
ied. Synchronization is referred to as a process wherein two (or many) dynamical 
systems are adjusted to a common behavior as time goes to infinity, due to coupling 
or forcing. Therefore, before discussing synchronization, a basic problem on continu- 
ation of the solution must be solved: For given initial conditions, can the solution of 
IT! ; 

Q \ coupled dynamical systems be extended to the infinite interval [0, +00)? In this paper, 

_C ■ we propose a general model of coupled dynamical systems, which includes previously 

studied systems as special cases, and prove that under the assumption of QUAD, the 
solution of the general model exists on [0, +00). 
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p : 1 Introduction 

o ■ 

In past years, collective behaviors of coupled dynamical systems have been widely stud- 
ied. In particular, synchronization in networks of coupled dynamical systems, as one of 
the simplest and most striking behaviors, has attracted increasing attention in mathemati- 
cal and physical literatures because of its potential applications in various fields, such as 
communication [1], seismology [2], and neural networks J3). 
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The word "synchronization" comes from a Greek word, which means "share time". To- 
day, in science and technology, it has come to be considered as "time coherence of different 
processes". Since the first observation of synchronization phenomenon was made by Huy- 
gens H| in the 17th century, many different types of synchronization phenomena have been 
found, e.g., phase synchronization, lag synchronization, full synchronization, partial syn- 
chronization, almost synchronization, and so on. In mathematics, synchronization can be 
defined as a process wherein two (or many) dynamical systems adjust a given property of 
their motion to a common behavior as time goes to infinity, due to coupling or forcing (see 
flU). For example, full synchronization requires that the difference between any two nodes 
converges to zero as time goes to infinity. Therefore, it is natural to raise following question: 
For given initial conditions, can the solution be extended to the infinite interval [0, +00)? 

For example, in the paper [1 1 J, the following coupled systems with a delay is considered: 



x\t) = f(x\t)) + c a v T 



x j (t-r) -x%t) 

2 = 1, ...,m, (1) 



where x' l (t) = [x\(t), ■ ■ ■ ,x t n (t)} T 6 M. n denotes the n-dimensional state variable of the 
i-th node, i — 1, . . . , m; f : W l — ► M. n is a differential function of the intrinsic system; c is 
the coupling strength; T = diag{7i, • • ■ , 7„} is the inner connection diagonal matrix with 
7i > 0, i = 1, . . . , n; > 0, for all i 7^ j, is the coupling coefficient from node j to node 
i; and r > is the coupling delay. It is assumed that J2f=i j^i a ij = 1> a a = ~ 1> f° r a ^ 
i = 1, . . . , m. And the following theorem was proved. 



Proposition 1 Suppose that there are a positive definite diagonal matrix P = diagjpi , • • ■ , p n } 
and a diagonal matrix D = diag{di, • ■ ■ , d n }, such that 

(x - y) T P[f(x) - f(y) -Dx + Dy] < -a{x - y) T (x - y) 

holds for some a > 0, any 1,1/6 R". Then, for sufficiently large coupling strength c and 
sufficiently small delay r, the coupled system ((7]) will be globally synchronized. 
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Here, a prerequisite condition in discussing synchronization is that the solution x l (t), 
i = 1, • • • , m, can be extended to the infinite interval [0, +00). However, in most papers on 
synchronization of coupled systems, such as E |71 [H dU and others, it is always assumed 
that for each initial condition, the coupled system under consideration has a unique solution 
for all time t > without any theoretical justification. 

In this short paper, we address this issue and propose a general model of coupled dy- 
namical systems, which includes previously studied systems as special cases. We prove that 
under the assumption of QUAD (Assumption (A5) in Section [2]), the solution of the general 
model exists on [0, +00). The assumption of QUAD is often used when using a Lyapunov 
function with a quadratic form to investigate the global synchronization (e.g., in Proposition 
1, and in iHHMTOl). Therefore, the theorem proved in this paper provides a theoretical basis 
for the discussion of synchronization of the coupled systems. 

The rest of the paper is organized as follows: In Section [2l we propose a general model 
of coupled dynamical systems. In Section [3l we present some fundamental theorems of 
retarded functional differential equations with infinite delay, which are taken from |fl4|. In 
Section |4l the main theorem is proved. We conclude the paper in Section [51 



2 Model descriptions 

In this section, we investigate the coupled dynamical systems described by the following 
retarded functional integro-differential equations: 

x\t) = f(t,x\t)) 

m poo 

+ y2 a ij( t ) / g(t, xJ (t - T ij(t) - s))dK i:j (s), 

3=1 J ° 

2=1,2, . ..,771, (2) 

where "'" represents the right-hand derivative, m is the network size, x l (t) E W 1 is the 
state variable of the 2-th node, t E [0, +00) is a continuous time, / : [0, +00) x W 1 — ► W 1 
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describes the dynamical behavior of each uncoupled system, A(t) = (a^(£)) G ]R mxm is 
the time- varying coupling matrix, which is determined by the topological structure of the 
network, g : [0, +00) x M™ -> f " is the output function, dKij(s) is a Lebesgue-Stieljies 
measure for each i, j = 1, . . . , m, and satisfies J °° |d-£Ty(s)| < +00. 

In addition, the following assumptions are necessary in discussion of retarded sys- 
tems: 

(Al) f(t, u) is continuous, and locally Lipschitz continuous with respect to u, i.e., in each 
compact subset W of [0, +00) x R", there exists a constant l(W) > such that 
\\f(t, u x ) - f(t, u 2 ) II < l(W) \\ Ul - u 2 1| for any (t, u k ) G W, k = 1, 2; 

(A2) A{t) = (a i:/ (t))™ =1 is continuous; 

(A3) g(t, u) is continuous, and there exists a continuous function nit) : [0, +00) — > IR + , 
such that \\g(t, Ui) — g(t, u 2 )\\ < K,(t)\\ui— u 2 \\ foranyt G [0, +00) and-U!,w 2 G M n ; 
(A4) For each z, j — 1, . . . , m, Tij(t) is continuous and nonnegative; 

(A5) There are a symmetric positive definite matrix P and a diagonal matrix A = diag{5i , . . . , 5 n } 
such that f(t, u) G QUAD(A, P), where QUAD(A, P) denotes a class of continu- 
ous functions h(t, u) : [0, +00) x IR n — > ]R n satisfying 

(«! - u 2 ) T P{[h(t, u x ) - h{t, u 2 )\ - A[ui - u 2 ]} 

< -e(ui - u 2 ) T (ui - u 2 ) (3) 

for some e > 0, all u\,u 2 G E n and t G [0, +00). 

Here, || ■ || can be any norm in R™ (Without loss of generality, in this paper we assume that 
|| • || is 2-norm). 

The model © includes many previously studied systems as special cases. In the follow- 
ing, we present several examples. 

Example 1 dKij(s) = 5(s), where 5(s) is the Dirac-delta function, i.e., 5(0) = 1 and 
5(s) = for s 7^ 0; A(t) = A is a constant matrix with zero-sum rows and nonnegative 
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off-diagonal elements; g(t,u) = Tu, where V is a constant matrix; Tij(t) = for each 
i, j — 1, . . . , m and all t > 0. Then, (|2|) reduces to the system with undelayed, constant and 
linear coupling discussed in /|6] UO\l : 

m 

x l (t) = f(t, x l {t)) + a^Tx^t), i = l,2,...,m. 

3=1 

Example 2 dKij(s) = 5(s), where S(s) is the Dirac -delta function; A{t) is a time-dependent 
matrix with zero-sum rows and nonnegative off-diagonal elements; g(t, u) = T(t)u, where 
V(t) is a time-dependent matrix; Tij(t) = for each i,j = 1, . . . , m and all t > 0. Then, 
(|2]) reduces to the system with undelayed, time-varying and linear coupling discussed in 

/Eliza/: 

rn 

% — 1,2, ... ,m. 

Example 3 dKij(s) = 5(s), where 5(s) is the Dirac-delta function; f(t,u) = f(u), i.e., 
f is independent oft; A(t) = A is a constant matrix with zero-sum rows and nonnegative 
off-diagonal elements, and satisfies an = —cfor i = 1, . . . , m; g(t, u) = Tu, where Y is a 
diagonal matrix with nonnegative diagonal elements; it) = r for i ^ j and Tu(t) = Ofor 
i = 1 , . . . , m. Then, ([2]) reduces to the system with delayed, constant and linear coupling 
discussed in /I771/: 

m 

x \t) = f{x\t)) + aijF[x j (t ~r)- x\t)l 

i = 1,2, ... ,m. 

Besides Examples 1-3, the model © includes coupled dynamical systems with nonlin- 
ear coupling, time-varyingly delayed coupling, distributedly delayed coupling, etc. 
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3 Preliminaries 



In this section, we present some fundamental results of retarded functional differential equa- 
tions with infinite delay, which will be used in the sequel. 

Firstly, we introduce some notations and definitions. 

Denote BC((—oo, a],R N ) the family of continuous functions <p mapping the interval 
(—00, a] into R N such that ||0|| =: sup{||0(0)|| : —00 < 9 < a} is finite. Also, denote 
C°°((-oo,a],R N ) = {(f) G BC((-oo,a),R N ) : lim«_ f _ 00 0(0) exists in R N }. When a = 0, 
we generally denote C°° = C°°((-oo,0],R N ). For a G R, B > 0, x G C°°((-oo, a + 
B], R N ), and t G [a, a + B] , we define x t G C°° as x t (9) = x(t + 9),6 G (-00, 0]. Assume 
is an open subset of R x C°°, h : Vt ^ R N is a given function, and " ' " represents the 
right-hand derivative; then, we call 

x{t) = h{t,x t ) (4) 
a retarded functional differential equation with infinite delay on VL. 

Definition 1 A function x is said to be a solution of Equation (I?]) on the interval I = 
[a, a + B) if there are a G R and B > such that x G C°°((-oo, a + B),R N ), (t, x t ) G Vt 
and x{t) satisfies Equation (0) for t G /. For given a G R, </? G C°°, z/a solution x of 
Equation (I?]) w defined on an interval [a, a + S), S > 0, and satisfies x a = if, then x 
is called a solution of Equation (I?]) with initial value ip at a or simply a solution through 
(a, <p). 

Definition 2 Suppose x(t) and y(t) are solutions with the same initial condition and satis- 
fies Equation (I?]) respectively on the intervals I and J whose left end points are 0. If I is 
properly contained in J and x{t) = y{t)for t G /, we say y is a continuation of x. If x has 
no continuation, it is called a noncontinuable solution, or a maximal solution. 
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Definition 3 We say h(t, 4>) is Lipschitz in <p in a compact subset W o/R x C°° if there a 
constant I > such that, for any (t, 4>h) £ W, k — 1, 2, 

\\h(t,cf )1 )-h(t, ( f> 2 )\\ <Z||0i-0 2 ||. (5) 

The following three lemmas on existence, uniqueness, and continuation of the solution 
of Equation (HJ), are used in the proof of the main theorem in the next section. The details 
can be found in lfl4l . 

Lemma 1 (Existence) Suppose Q is an open subset in R x C°° and h : Q —>■ ~R N is contin- 
uous. Then, for any (a, (p) G Q, there exists a solution of Equation (0) through (a, (p). 

Lemma 2 (Uniqueness) Suppose Q is an open subset in R x C°° and h(t, <fi) is Lipschitz 
in (p in each compact subset of Q. Then, for any (a, ip) G Q, there exists at most one 
noncontinuable solution of Equation 0j) through (a, tp). 

Lemma 3 (Continuation) Suppose Q is an open subset in R x C°°, h : Q — > M. N is con- 
tinuous, and x is a noncontinuable solution of Equation (@) defined on I = [a, a + B). 
Then, for every compact subset W ofVt, there is a ty/ in I such that (t, x t ) ^ W for all 
t G (t w ,a + B). 



4 Main result 

In this section, we prove the following theorem. 

Theorem 1 Suppose that Assumptions (Al )-(A5) hold. Then, for any (p(6) = [ip 1 (8) T , . . . , ip m (6) T ] T 
withtp l (0) G oo, 0], R"), there is a unique noncontinuable solution x(t) = [x 1 ^)' 7 , . . . ,x m (t) T ] 

of Equation (HI) through (0,ip). Moreover, the interval of existence of the solution x is 
[0,+oo). 
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Proof : By Assumptions (A1)-(A4) and Lemmas 1-2, it is clear that for the integro- 
diffential system ©, there exists a unique noncontinuable solution x(t). In the following, 
we will prove that the interval of existence of the solution x(t) is [0, +00). 

We employ "proof by contradiction", and suppose that the interval of existence of the 
noncontinuable solution x(t) is [0, b), where b is a positive constant. 

Firstly, by Assumptions (A1)-(A4), we can find positive constants a, (3 and 7 such that 

Hfl^jUi) -g{t,u 2 )\\ < a\\ui -u 2 \\ 
holds for all m, u 2 E R n and t E [0,b), and 
ki(*)| <P, 

m /»oo 

||/(*,x*(o)) + J2 a ij(t)g(t,x j (a)) / d%( s )|| < 7 . 
3=1 Jo 

hold for all i,j — l,...,m and t E [0, b). 

Now, we will show how the assumption of QUAD (Assumption (A5)) plays an important 
role in the proof. 

Since f(t, u) E QUAD(A, P) (Assumption (A5)), it is clear that there is a constant 
5 > such that for all U!,u 2 e W 1 and t > 0, 

(«i - u 2 ) T P[f{t, Mi) - f(t, u 2 )\ <8{u x - m 2 ) T (mi - u 2 ). 

Denote 

25 + 2a3\\P\\K 2m^\\P\\ 

>0, 



P 
min 



where m is the number of the nodes, ||P|| is the 2-norm of the matrix P, A^ in is the minimum 
eigenvalue of the matrix P, and K = YlT=i Sj=i Io° \dKij(s)\- 

Since the matrix P is symmetric positive definite, we can define a norm in M" m : 



|x(*)||p = (^x i (t) T Px i (t) 
1=1 
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and two nonnegative functions: 



V(t) = ^\\x(t)-x(0)\\ 2 P , 




x(s)-x(0)\\ f 



p > 



Clearly, V{t) < M (t).We claim that Af(t) < Af (0)e* for all t G [0, 6). 

In fact, at any t G [0, 6), there are two possible cases: 

Case 1: V(to) < M(t ). In this case, by the continuity of \\x(t) — x(0)\\ 2 P , M(t) is non- 
increasing at to. 

Case 2: V(t ) = M(t ). 

Calculating the right-hand derivative of V with respect to time along the trajectories of 
©, one has 



hi 



V(t ) = ]T(x l (t ) - x\0)) J pj(to) 



i=l 



rn 



5>'(tb) - Ao)) t p f(t ,x%)) 



i=l 




HI 




i=l 





fit 



< 



5j2(x\to) ~ x l (0)) T (x l (t ) - x l (Q)) 



i=i 
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+EEM*o)|||p||||^o)-^(o)|| 



i=i i=i 

"OO 



/•OO 

x / \\g(to, X 3 (t - Tij(t ) - s)) 

JO 

-^o,^(0))|||d^-( S )| 

m 

+Y.\\ p \\¥^~ xl ^Mf^ xl ^)) 



i=l 
m 



ItL /»oo 

+ y^a ij {t Q )g{t G ,x j {Q))) / di^(s 

m m 

5\\x(t ) - x(0)\\ 2 + ^||P||EEH'°) " ^(0)|| 

i=i j=i 



X 



oo 



||^(t -^(to)-s)-^(0)|||d^( S )| 

m 

+7WIEIW*o)-^(0)|| 



1=1 

|2 



5||x(t ) - x(0)|| 2 + a(3\\P\\ EE N*°) " 

i=i i=i 



oc 



|ar(*o - 7ij(* ) -s) -x(0)|||d^(s)| 



x 

+m7||P||||x(i ) — ar(0)|| 
*\\ x (to) - x(0)\\ 2 p 

min 

/ill pll m m 



X 



/•oo 

/ ||x(t - Tijito) -s)- x(0)\\ p \dK i:j {s)\ 
Jo 



+ ! ^ / =N*o)-*(0)|| p -l 
V \ 

5 



p 

min 



-f 



2M(t ) 



^v^v^EE / l d ^-( s )l 

min i=1 j=1 JO 
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25 + 2am\P\\K 2m 7 Pi , 
J > " " + -7= W*,, 

min V A min 

^(to) 



In summary, we conclude that M(t) < M(0)e vt for all t E [0, b), which implies V(t) < 
M(0)e^ and 

Mt)\\ < -^=\\x(t)\\p 

V A min 

< ^=(||x(0)|| P +||x(t)-x(0)||p) 

V min 

= -7^(IKo)||p + V2y(t)) 

V A min 

< -^(||x(0)|| P + V2M(0)^) 

V A min 

< -yL=(\\x(0)\\p + \/2M(0)e<t>) (6) 

V min 



for all t E [0,6). 

Now, pick a compact set 

W= S.(t,ip) GRx C°°((-oo,0],E nw ) 

< max 



< t < b, and 



-7^(ll^(0)||p + V2M(0)e^), |M 

V A n 



min 



where (/? is the initial value. By the inequality ©, we conclude that (£, x t ) E W for all 
t E [0,b), which contradicts Lemma 3. 

Therefore, the interval of existence of the noncontinuable solution x is [0, +00). Theo- 
rem is proved completely. 
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5 Conclusions 



In this paper, we propose a general model of coupled dynamical systems, which includes 
previously studied systems as special cases, and prove that under the assumption of QUAD, 
the solution of the general model exists on [0, +00). 
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